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Abstract. 

Although it is weU known that the solar acoustic mode frequency increases as the solar 
activity increases, the mechanism behind it is still unknown. Mode frequencies with 20 < I < 900 
obtained by applying spherical harmonic decomposition to MDI full-disk observations were used. 
First, the dependence of solar acoustic mode frequency with solar activity was examined and 
evidence of a quadratic relation was found indicating a saturation effect at high solar activity. 
Then, the frequency dependence of frequency differences between the activity minimum and 
maximum was analyzed. The frequency shift scaled by the normalized mode inertia follows a 
simple power law where the exponent for the p modes decreases by 37% for modes with frequency 
larger than 2.5 mlfz. 



1. Introduction and Data 

It seems that the responsible mechanism for the mode frequencies change with solar activity is 
restricted to the outer layers of the Sun. However, at the moment, there is no general agreement 
in the precise physical cause. The frequency variations of medium and high-/ modes with solar 
cycle were analyzed to help determine its physical origin. The solar-radio 10.7-cm daily flux 
(NGDC/NOAA) was used as the solar activity proxy. 

The mode frequencies were obtained by applying spherical harmonic decomposition to MDI 
full-disk observations during the Dynamics and Structure observing modes. The first one has 
higher spatial resolution and is available every year for two or three months of continuous data. 
Data for 1999 to 2008 were used. The second one consists of 72-day time series from early 1996 
to April 2010 (Fig. [1]). They will be called the Dynamics and Structure sets from now on. 

The mode frequencies were obtained using two distinct methods. The medium-/ method ([5] 
and [7j) and the high-/ method ([B]) which is applied only to the Dynamics set (Fig. [T] left panel). 
The second method is used when the spatial leaks of the modes overlap with the target mode 
making it more difficult to estimate the mode frequency. 



2. Mode frequency variation vi^ith the solar activity cycle 

There is a very high linear correlation of the mode frequency variation with several solar activity 
indices. However, deviations from a simple linear relation has been reported (see, for example, 
[2]). In this analysis, there is some indication of a quadratic relation: 



fg(n,/, Fr) = co{n,l) + ci{n,l) ■ Fr + C2(n, /) • F^ 



(1) 
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Figure 1. Right panel: Mode coverage. Modes obtained by the medium-/ method are in red 
and by the high-/ method in black. Left panel: Mean radio flux for the observational periods in 
the Dynamics set (in black) and in the Structure set (in red). 



where Fr is the relative activity index defined as the mean solar-radio flux of a given observational 
period divided by the maximum mean flux of the set of observations (Fig. [2|). The quadratic 
polynomial seems to indicate a saturation effect at high solar activity. A similar effect has been 
seen in frequencies at activity regions with a large surface magnetic field using ring analysis 
([!]). The saturation occurs when ci + 2 ■ C2 ■ = 0, i.e.: = — ci/(2 • C2). 
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Figure 2. Examples of frequency variation with the relative activity index for different modes. 
The top left panel is using the Structure set and the others using the Dynamics set. 

Fig. [3] shows the coefficient C2(n,/) as a function of ci(n, /). For the Structure set, the slope 
of the C2 — ci linear relation is —0.3077 it 0.0025 with a zero y-intercept (blue line in both panels 
of Fig. [2D, indicating a saturation at (357.2 it 2.9) x 10~^^ J s~^ m~^ Hz~^. For the Dynamics 



set, the slope is -0.2813 ± 0.0083 with a very small y-intercept (-0.0169 ± 0.0052 /iHz). For 
comparison, the linear regression was performed on the same ci range as the Structure set: 
< ci < 4.34 /iHz. Modes obtained with both medium- and high-/ methods were included in 
the fitting. Only 12% (133) of the modes in the Dynamics set were obtained with the high-Z 
method. The saturation occurs at (329.2 it 9.7) x 10^^^ J s~^ m~^ Hz~^. This is more or less in 
agreement (differing by 2. 80") with the value estimated by the Structure set. The slope increases 
as ci (and |c2|) increases (Fig. [H]). For ci > 4.34 |uHz, the slope is —0.455 ± 0.028 and the 
y-intercept 0.50 ± 0.29 /xHz, giving a saturation at (204 ± 12) x lO'^^ j g-i ^^-2 Hz^i. 




-20 ^ i - 

-20 -10 10 20 30 
Ci (/xHz) 




1 2 3 4 5 
Cj (/iHz) 



Figure 3. Left panel: Dynamics and Structure sets. Right panel: Structure set. The green 
circles in the left panel are the same as the black circles in the right panel. The red and 
black circles in the left panel are for the Dynamics set using the medium- and high-/ methods 
respectively. The light-blue line is the mean slope of the simulations for the Dynamics set. 

To test the significance of the quadratic relation in the data, the observed frequencies were 
fitted using a linear relation, iyi{n,l,Fr) = do{n,l) + di{n,l) ■ Fr, and noise was added to the 
fitted frequencies. To each fitted frequency i'i{n,l, Fr), the noise was generated as normally 
distributed random numbers with the same standard deviation as the error in the observed 
frequencies. For each (n, /) mode and each realization, a quadratic polynomial was fitted. Fig. H] 
shows the coefficients obtained for one realization and the ones obtained for the Structure 
set for comparison. A straight line was fitted to C2(ci) for each realization over the entire 
range of ci of each set. The slope averaged over 1000 realizations is —0.01437 zt 0.00011 
and —0.712222 it 0.000060 for the Structure and Dynamics set respectively. The first one is 
represented by the black line in Fig. H] and the second one by the light-blue line in the left 
panel of Fig. [31 These slopes are very different from the ones obtained by the corresponding set 
indicating that the quadratic relation is not due to noise but to a real signal. 

For the Dynamics set, there are several modes with a negative coefficient ci (Fig. [3|). Nine 
percent (195) of the modes have ci < — e^, where is the fitting uncertainty. For the Structure 
set, there are two percent (27) of the modes with ci < — ec^. An example of the frequency 
variation with solar cycle for one of these modes is in the bottom-right panel of Fig. [2j The 
frequency of this mode is approximately constant for Fr < 0.8. To see if this is true for all the 
modes with a negative ci, the mean relative frequency at a given time t is defined as: 

//? (f\\ = J- V ^n.;M - I^nA'^M^r)] . . 
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Figure 4. Coefficient C2 versus ci for 
one realization (red circles) and from the 
observations using the Structure set (black). 
The black circles are the same as the black 
circles in the right panel of Fig. [31 The 
black and blue lines are the linear fit to the 
simulation and to observations respectively. 
The blue line is the same as the blue line in 
Fig.H 



where is the number of modes. (i?^(t)) was estimated over all modes with ci > and 
ci < —Cc^ given by the black and red circles in Fig. [5] respectively. For the Structure set, 
N = 1573 for ci > and N = 27 for ci < — 6^. For the Dynamics set, N = 1272 for ci > 
and = 60 for ci < — Cc^. Only modes with a Spearman's correlation coefficient larger than 
0.75 were used in the average of the Dynamics set to decrease the scatter. Modes with a negative 
ci seem to be less sensitive to the solar cycle, more so away from the solar maximum. Most 
of these modes have a small I'/L, between 3 and 25 /iHz. In the Dynamics set, most of them 
that are p modes have 3 < u < 3.6 mHz and / > 200 while the / modes have / > 600. In the 
Structure set, most that are p modes have u < 1.6 mHz and / < 170 while the / modes have 
100 < I < 140. It is not clear why these modes are less affected by the solar cycle. If this 
is proved to be true and not an artifact of the data, it will help in the understanding of the 
mechanism or mechanisms responsible for the mode frequency variation with solar activity. 
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Figure 5. Left panel: Structure set. Right panel: Dynamics set. Modes with ci > are in 
black and modes with ci < —e^ in red. The relative activity index is in blue. In the right panel, 
the medium- and high-Z modes are given by the dashed and dotted lines respectively. 

From Fig. O (left panel), the minimum between solar cycle 23 and 24 happened around August 
2008. The frequencies were at a minimum from November 2007 to August 2009. As pointed 
out by other authors (see [4] and references within), the mode frequencies were smaller than the 
previous minimum. The frequency averaged over all medium-Z modes is ~ 10% of the minimum- 
to- maximum variation smaller than in the previous minimum: (i?j^(July 1996)) — (i?;^ (August 
2008)) ^ 0.1. 



3. The minimum-to-maximum frequency shift 

Here, the minimum-to-maximum frequency shift 5i'{n, I) is defined as the difference between the 
fitted frequency given by Eq. (1) at the maximum and minimum solar activity. The minimum 
and maximum activity was chosen to be those of the Dynamics set, which are around April 
2008 and May 2002 respectively (Fig. [1] right panel). The differences in 6v{n,l) using a linear 
regression instead of Eq. (1) are small. The standard deviation of the differences is 0.0031 ;uHz 
and 0.25 ;uHz for the Structure and Dynamics set respectively. Most of the difference is due to 
the high-/ modes in the Dynamics set. The effect of these differences in the estimated a and 7 
coefficients described below is in their 4-th significative digit or higher. 

Most of, if not all, the perturbation to the structure of the Sun associated with the solar 
cycle is believed to be located in a thin layer near the solar surface which leads to a frequency 
perturbation roughly proportional to the inverse of the mode inertia InX- ^^n,i oc v^JIni- The 
value of a depends upon the physical mechanism responsible for affecting the mode frequencies 
during the solar cycle. Some of the predicted values are -1, 1 and 3 (see [3] and references 
within). Fig. [6] (right panel) shows 8v{n, I) scaled by for both sets. Only modes with ci > 
were plotted. Fitting only the p modes, a varies from approximately -4 for u < 1.6 mHz, to 0. 
for 1.66 < 1/ < 2.5 mHz, to 1.43772 ± 0.00018 for 2A < u < 3.55 mHz, to 6.417 ±0.023 for 1/ > 4 
mHz. The variation of a indicates a more complex mechanism than anticipated. 




Figure 6. Only modes with ci > are plotted. 

Recently, [6] observed that scaling the frequency shift with the mode inertia normalized 
by the inertia of a radial mode of the same frequency {Qn,i) also follows a simple power law: 
Sun^i (X i^^i/Qn,ii but with a constant exponent. Extending this previous analysis to include 
data from the end of solar cycle 23 gives similar results (Fig. [6] left panel). However, there is 
a clear change in the coefficient 7 around 2.5 mHz. For the p modes, 7^ = 6.04466 it 0.00030 
for I' < 2.5 mHz and 7p = 3.78657 it 0.00019 for > 2.5 mHz. In the previous analysis, it was 
estimated 7^ = 3.60 it 0.01 at all frequencies and it was argued that around 2.3 mHz there was 
a step in the p-mode frequency. This analysis show clearly that there in fact a change in 7^. 

The p modes obtained by the high-/ method show a slightly different slope than those obtained 
by the medium-/ method. In Fig. [6l the p modes obtained by the high-/ method are concentrated 



at log[z/(//Hz)] > 3.55 where they account for 70% of the modes (including the low n modes 
indicated by the colored stars). Although an artifact from the analysis can not be discarded 
due to the difficulty in estimating unbiased mode frequencies at high-/ and/or high-frequency, 
it could be an indication of a degree dependence in the scaled frequency shift. 

It is expected that the / modes are affected by the solar cycle in a different way than 
the p modes since they have very different properties. Fig. [7] shows the /-mode frequency 
shift. The high-Z / modes in the Dynamics set has a distinctive behavior (red stars in Fig. O 
with 3.25 < log[i^(/uHz)] < 3.4). Fitting only the modes obtained by the medium-/ method, 
jf = 7.0880 lb 0.0019 which is 17% larger than 7^ at the corresponding frequency range. 




Figure 7. Only f modes. In red are the Dynamics set and in black the Structure set. The 
high-/ modes are those with 1^ > 1760 fjMz. The blue line is for 7^ = 7.0880 it 0.0019. 



4. Summary 

There is some evidence of a quadratic relation of the frequency shifts with the solar-radio flux 
indicating a saturation effect at high solar activity. The saturation occurs at ~ 350 x 10^^^ J 
s~^ m~^ Hz~^ and decreases as ci increases. 

The minimum-to-maximum frequency shift scaled by the mode inertia is proportional to z/° 
where ap increases from -4 at low frequency to 6.4 at high frequency for the p modes. Scaling 
by the normalized mode inertia, there is a sharp change in the p-mode exponent at 2.5 mHz. 

Acknowledgments 

The frequencies obtained by the medium-/ method were calculated by J. Schou and T. Larson. 
References 

[1] Basu, S., Antia, H. M. and Bogart, R. S. 2004, ApJ, 610, 1157 

[2] Chaplin, W.J., Elsworth, Y., Miller, B.A. and Verner, G.A. 2007, ApJ, 659, 1749 

[3] Gough, D.O. 1990, Lecture Notes in Physics, 367, 283 

[4] Jain, K. et al. 2010, ASP Conf. Ser, 428, 57 

[5] Larson, T. and Schou, J. 2009, ASP Conf. Ser., 416, 311 

[6] Rabello-Soares, M. C, Korzennik, S. G. & Schou, J. 2008, Solar Phys., 251, 197 
[7] Schou, J. 1999, Ap,J, 523, L181 



